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| The history of linear differential equations is over 350 years. By using Frobenius 

method and putting the power series expansion into linear differential equations, 
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the recursive relation of coefficients starts to appear. There can be between two 



and infinity number of coefficients in the recurrence relation in the power series 
expansion. During this period mathematicians developed analytic solutions of only 
two term recursion relation in closed forms. Currently the analytic solution of three 
term recurrence relation is unknown. In this paper I will generalize the three term 
recurrence relation in the linear differential equation. This paper is 2nd out of 10 
in series "Special functions and three term recurrence formula (3TRF)". The next 
paper in series deals with the power series expansion in closed forms of Heun function 

m ■ 

by Choun |arXiv:13 03.0830j. The rest of the papers in the series show how to solve 
mathematical equations having three term recursion relations and go on producing 

X 

the exact solutions of some of the well known special functions including: Mathieu, 



Heun, Biconfluent Heun and Lame equations. See section VII for all the papers and 
short descriptions in the series. 
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I. INTRODUCTION 



Mathieu functions 6 , is an example of three term recurrence relation appears in physical 
problems involving elliptical shapes 7 or periodic potentials, were introduced by Mathieu 
(1868)- when he investigated the vibrating elliptical drumhead. 

Mathieu function has been described using numerical approximations (Whittaker 1914 3 , 
Frenkel and Portugal 2001 1 ). Whittaker tried to obtain the analytic solution of Mathieu 
equation using Floquet's theorem and he reached the conclusion using three term recursive 
relation. He did not represent the solution of Mathieu equation in closed forms because 
of its three term recursive coefficients. He argued that "While the general character of the 
solution from the function theory point of view is thus known, its [Mathieu] actual analytical 
determination presents great difficulties. The chief impediment is that the constant \i (the 
Mathieu exponent) cannot readily be found in terms of a and q."—. a and q (see (1) in Ref. O) 



are corresponding to A and — |g (see (1) in Ref. Il5l ). In my opinion, Mathieu functions are 
difficult to represent in analytic closed forms and in its integral forms because of the three 
recursive coefficients. 

Heun function, is an example of three term recurrence relation, generalizes all well-known 
special functions such as: Spheroidal Wave, Lame, Mathieu, and hypergeometric 2 -Pi, 
and qFi functions. The Heun equation has four kinds of confluent forms such as Confluent 
Heun, Doubly-Confluent Heun, Biconfluent Heun and Triconfluent Heun equations. Mathieu 
equation is the special case of Confluent Heun equation. 

According to Hortacsu, 'Heun differential equation^ starts to appear in astronomy or 
in general relativity 10 inevitably. Heun equation is a general equation of all well-known 
special functions; Mathieu, Lame and Coulomb spheroidal equations. The power series 
expansion on Heun function can not be described as two term recursion relation any more. 
The coefficients in a power series expansions of Heun equation have a relation between three 
different coefficients. - 

For the past 350 years, we have been only using in two term recurrence relation. Due 
to its complexity more than three term case has been neglected. However, "since 1930, 
we do not have simple problems to solve in theoretical particle physics and scientists and 
mathematicians doing research on this field have to tackle more difficult problems, either 
with more difficult metrics or in higher dimensions. Most of the difficult problems must 



include three term or more" (Hortacsu 201 1 5 ). 

In this paper I will construct three term recurrence relation in the form of power series 
expansion for the cases of polynomial and infinite series. In the following papers (see section 
VII) using the method described in this paper I will show how to obtain exact solution (1) in 
the power series expansion, (2) in the integral formalism and (3) in the generating function 
of any linear ordinary differential equations having three term recurrence relation. Using 
the method described in this paper one might be able to obtain exact solutions of higher 
term recurrence formula; four, five, • • • , m , oo term recurrence formulas. If this is possible, 
then one would be able to generalize all homogeneous and inhomogeneous linear differential 
equations. 

Most problems in nature turns out to be nonlinear. For simplification purposes we usually 
try to linearize these nonlinear system. The systems are linearized using certain methods 
of simplification resulting in better approximation of physical systems. All physics theo- 
ries (E&M, Newtonian mechanics, quantum mechanic, QCD, supersymmetric field theories, 
string theories, general relativity, etc), generally involve solutions of linear differential equa- 
tions. Unfortunately in some instances there are no analytic solutions to these physical 
systems. Hopefully the methods and procedures described in these papers (see section VII) 
will help obtain better analytic solutions to linearized nonlinear systems. 



II. LAME EQUATION, FROBENIUS METHOD AND THREE TERM 
RECURRENCE RELATION 

There are no exact solutions in closed forms for second order (or higher order) ordinary 
linear differential equations consisting of three term recurrence relation: some of the exam- 
ples are Lame function, the generalized Lame function, Heun's equation, G.C.H. function^, 
Mathieu function, etc. Two term recurrence relation are solvable. Some examples are: Leg- 
endre function, hypergeometric function, Kummer function, Bessel function, etc. Let's think 
about an example which has no analytic solution in closed forms in detail. 

&y 1 f 1 1 J_\ dy _ a(a + l)t + /3 

dt 2 2 \t-a t-b t-c) dt t(t - a)(t - b)(t - cf U 

([1]) is Lame differential equation. If a is not positive integer, the solution of it is called as 
the generalized Lame function. Replace t by x + a in ([!]) . By using the function y(x) as 
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Frobenious series in it, 

CO 

y(x) = c n(t ~ «) n+A = E C " xn+X ( 2 ) 

n=0 n=0 

Plug in (PQ). And its recurrence relation is 

c n+ i = j4„ c n + S n c„_i ; n > 1 (3) 



where, 



K n = A n + — ; n > 1 (4a) 

A, 



n-1 



= — = — (4b) 

{a(a + l)a + ^} - 2 2 (2a - 6 - c)(n + A) 2 
n ~ 2 2 (a - 6) (a - c)(n + A + l)(n + A + ±) ( C) 

{ a + 2(r2 + A)-l}{ a -2(r2 + A) + 2} 
- 2 2 (a-6)(a-c)(n + A + l)(n + A + |) 1 J 

All other differential equations having no analytic solution in closed forms can be described 
as in ([3]). If I get a formula of ([3]) type, I can apply it to all other functions having no 
analytic solution in closed forms to diverse areas such as Lame function, generalized Lame 
function, Mathieu function, Heun function, G.C.H. function^, etc. 



III. INFINITE SERIES 

Assume that 

ci = Aq c (5) 

(jSJ) is a necessary boundary condition. The three term recurrence relation in all differential 
equations having no analytic solution in closed forms follow (JSJ). 

<jj— i 

IT ^ n = ^ where a; is positive integer including (6) 

7l=d; 

([6]) is also a necessary condition. Every differential equations having no analytic solution in 
closed forms also take satisfied with (jSJ). 
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My definition of refer to BiBjBkBi. Also, Aj,fc,« refer to AiAjA^Ai. For n = 

0,1,2,3,-.., © gives 

c 

Cl = A c 

c 2 = (A ,i + B l )c 

c 3 = (A,i,2 + A 2 B X + A 5 2 )co 

C4 = (A),l,2,3 + ^2,3^1 + A ,3-B 2 + A ,i£ 3 + S li3 )co 

C5 = (^4.0,1,2,3,4 + A,3,4-8l + A),3,4-8 2 + Aq^Bs + ^9,1,2 -84 + A 4 S lj3 

+A 2 B 1A + A B 2A )c 

c 6 — (^4-0,1,2,3,4,5 + ^2,3,4,5-81 + ^0,3,4,5-82 + A), 1,4,5-83 + ^0,1,2,5-84 + ^-0,1,2,3-85 + ^4,5-81,3 

+^-2,5-8l,4 + ^-0,5-82,4 + ^-2,3-8l,5 + ^-0,3-82,5 + A), 1-83,5 + -8l,3,5)co 
C7 = (A,l,2,3,4,5,6 + ^2,3,4,5,6-8l + ^0,3,4,5,6-82 + ^0,1,4,5,6-83 + A), 1,2,5,6-84 + A, 1,2,3,6-85 

+ A, 1,2,3,4-86 + A, 5,6-81,4 + ^-4,5,6-8l,3 + ^0,5,6-82,4 + ^-2,3,6-8l,5 + ^0,3,6-82,5 

+ A, 1,6-83,5 + A,3,4-8l,6 + ^-0,3,4-82,6 + A), 1,4-83,6 + A,l,2-84,6 

+A-8i,3,6 + A 6 5i i3) 5 + A-81,4,6 + A S2,4,6)C0 
C8 = (A,l,2,3,4,5,6,7 + A,3,4,5,6,7-8l + A,3,4,5,6,7-82 + A,l,4,5,6,7-83 + A, 1,2,5,6, 7-84 + A,l,2,3,6,7-85 

+A, 1,2,3,4,7-86 + A, 1,2,3,4,5-87 + ^2,5,6,7-81,4 + A,5,6,7-82,4 + A,3,6,7-Bl,5 + A,5,6,7-8l,3 

+A,3,6, 7-82,5 + A,l,6,7-83,5 + ^2,3,4,7-81,6 + A),3,4,7-82,6 + A,l,4,7-83,6 + A,l,2, 7-84,6 

+A,3,4,5-8l,7 + A,l,4,5-83, 7 + ^40,3,4,5-82, 7 + A,l,2,5-84,7 + ^0,1,2,3-85,7 

+A,7-8l,4,6 + ^0,7-82,4,6 + ^4,5-81,3,7 + A,5-Bl,4,7 + ^44,7-8l,3,6 + ^0,5-82,4,7 

+A-6, 7-8l,3,5 + ^2,3-81,5,7 + A,3-82,5,7 + A,l -B 3j5 ,7 + -Bl,3,5,7)co 

(7) 

In (j7j) the number of individual sequence c n follows Fibonacci sequence: 1,1,2,3,5,8,13,21,34,55,- ■ ■ . 
The sequence c n consists of combinations A n and B n in (J7J). First observe the term inside 
parentheses of sequence c n which does not include any A^s: c n with even index (c , c 2 , 
c 4) - ■ •). 
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(a) Zero term of Ai's 



Co 

c 2 = Si c 
C4 = Bi i3 c 

(8) 

c 6 — -Bl,3,5 c 
C8 = -5l,3,5,7 c 



When a function analytic at x — 0, is expanded in a power series, we write 

CO CO 

y( x ) = C « x?1+A = $Z = y ° ^ + + V2 ^ + ' ( 9 ) 

n=0 m=0 

where 

CO 

y m (x) = ^cf^ +A (10) 

1=0 

A is the indicial root. y m (x) is sub-power series that have sequence c n including m term of 
Ai's in ([7]). For example yo(x) has sequences c„ including zero term of A n 's in ([7]), has 
sequences c n including one term of A n 's in ([7]), y 2 (^) has sequences c n including two term of 
A n 's in ([7D, etc. Substitute P) in ([10]) putting m = 0. 

co f n— 1 

y (x) = c J2\l[B 2il+1 L 2 " +A (11) 

n=0 Ui=0 J 



Observe the terms inside parentheses of sequence c n which include one term of A„'s in 
([7J): c n with odd index (ci, c 3 , c 5 ,- • ■ ). 

(b) One term of A n 's 
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Ci = A C 

C3={A (f)l+A 2 (f)fi 1 } Co 

c 5 = { A 1 + A 2 (*f ) Sl + A 4 (g£) B1i8 } Co 

c 7 ={A 1 + A 2 Bl + A 4 S 1i3 + A 6 (fg) 5,3,5)00 

c 9 ={A 1 + A 2 (^f) 2* + A 4 B 1>3 + A 6 (^m) B lj3 , s (12) 

en = { A Q 1 + A 2 (^f^) 5x + A 4 (^ff^) Si,3 + A 6 (^gf ) 5 M,5 

+As jBl ' 3 ' 5 ' 7 + Al ° (Sjejio) ^-3,5,7,9} CO 



( |T2l) is simply 



n f ii— 1 n— 1 ^] 

C2n+1 = I A 2il Yl B 2i 2 +1 JJ #2i 3 +2 > C (13) 

11=0 I. 12=0 «3=U J 



Substitute ( fTBl in ( |T0|) putting m = 1. 



yi(^) = co j2 j e ^ n ^ 2+ i n ^3+2 x 2 - +i+a (14) 

n=o Ln=o L «2=o «3=n J J 



Observe the terms inside parentheses of sequence c n which include two terms of A n 's in (JTj): 
c n with even index (c 2 , c 4 , c 6 ,- ■ ■ ). 
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(c) Two terms of A n J s 



C2 — A),1 C 



« = [a» {M (i) (§f) 1 + * (ft) (fe) 1} + * {* (ft) (fe) Br) }<:» 
A» (1) (%f ) 1 + A 3 (f ) 1 + ^ (¥) (f* 



{* (ft) * + * (%0 (is) M + * {* (fe) (is) M }<* 

Cs = U{a, (i) (ag*i) 1 + ,4 3 (f ) (%?f ) 1 + ^5 (^) (%5f ) 1 + a, (fgj 



+^{^ (ft) (%Jf ) B: + A ) (^) Bl + A, (Igjl) B, 



3.5 



(15) 



({15"]) is simply 



n ( n f ii — 1 12 — 1 re— 1 

C 2 n+2 = J]] ] ^2u \ ^2i 2 +l JJ S 2i3+1 JJ S ^+ 2 II B ^+ 3 }" }" C ( 16 ) 

11=0 L 12=11 L «3=0 14=11 i5=«2 



Substitute ( 1161) in ( flQl) putting m = 2 



oo ( re ( n ( ii — 1 12 — 1 re— 1 



h,(z) = c ^ \ a ^ e A 2i 2+ i n ^3 + i n n ^ 5+3 \ \ \ ^ +2+a a?) 

n=0 lu=0 L i2=h I «3=0 U=h *5=*2 



Observe the terms inside parentheses of sequence c n which include three terms of A n 's in 
(J7J): c n with odd index (c 3 , c 5 , c 7 ,- ■ ■ ). 
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(d) Three terms of A„s 
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(18) 



(|T8|) is simply 

re f n ( n ( ii — 1 22 — 1 13 — 1 re— 1 

C2n+3 = E ) ^ 2i l E ) ^2i 2 +l E S ^2i 3 +2 -S 2 j 4+ i B 2 i 5+ 2 Y\_ ^6+3 ^2i 7 +4 f ? ? r " 

il=0 I. «2=il I «3=«2 I J4=0 25=21 26=22 27=23 

(19) 

Substitute ( !T9|) in ( ITOj) putting m = 3. 

oo ( n f n ( n ( 22 — 1 

2/s(a;) = c ^ < E ] A 2il E < A 2i2+ i E i A 2 i 3+ 2Y[ B ^ +1 II B **+ 2 

n=0 V U=0 L «2=«i I «3=*2 I 24=0 25=11 



13— 1 n— 1 

X 

26=22 27=23 



n s 2J6+3 n B 2i7+i (20) 



By repeating this process for all higher terms of A's, we obtain every y m {x) terms where 
m > 4. Substitute (ITT]) . ([HI) . (FiTl) . (I2"0"j) and including all y m {x) terms where m > 4 into ([HD- 
Then general expression of y(x) for infinite series is 

oo 

V( x ) = ^c n x n+A = y (x) + j/i (x) + y 2 (x) + y-iix) H 

n=0 

{oo / re— 1 \ oo ( n ( re— 1 

e n - 2 - +A + e e ^ n n ^ ^ ^ ^ m+a 

n=0 \ii=0 / n=0 Ui=0 I i 2 =0 23=21 

oo ( oo f AT— 1 / n 

+E E II E 

JV=2 L n=0 L fe=l \i h =ik-l 

n ( JV-1 / «i+l-l \ n-1 ^ ^ ^ > 

X E ] A ^N + {N-l) II ( II B 2il+1+N+{ i +l) J JJ J B 2i2JV+1+(A r +1) ^ \ \x 2n+N+X \ 

2jV=2jV-l I 1=0 \il + l + N=k / *2JV+l=«iV ) ) ) ) 

where, i = (21) 



IV. POLYNOMIAL WHICH MAKES B n TERM TERMINATED 

Now let's investigate the polynomial case of y(x). Assume that B n is terminated at 
certain value of n. Then each yt(x) where i = 0, 1, 2, • • ■ will be polynomial. Examples of 
these are Heun's equation, G.C.H. function^, Lame function, etc. First B 2 k+i is terminated 
at certain value of k. I choose eigenvalue j3o which B 2 k+\ is terminated where /?o = 0, 1, 2, • • • . 
B 2 k+2 is terminated at certain value of k. I choose eigenvalue f3i which B 2 k+ 2 is terminated 
where j3% = 0, 1, 2, • • • . Also B 2 /, +3 is terminated at certain value of k. I choose eigenvalue 
(3 2 which B 2k+3 is terminated where j3 2 = 0, 1, 2, • ■ ■ . By repeating this process I obtain 

B Wl+{i+1) = where i, A = 0, 1, 2, • • • (22) 

10 



In general, the two term recurrence relation for polynomial has only one eigenvalue: for ex- 
ample, the Laguerre function, confluent hypergeometric function, Legendre function, etc. 
But three term recurrence relation has infinite eigenvalues which are fa, where i, fa = 
0,1,2,-... 

First observe the term in sequence c n which does not include any A n 's in (jSj): c n with 
even index (cq, C2, c±,- ■ ■ ). 



(a) As A)=0, then ^=0 in 



Co 



(23) 



(b) As po=l, then B 3 =0 in ©. 



c 

C 2 = 5iC 



(24) 



(c) As /3 =2, then B 5 =0 in 



Co 

c 2 = Sico (25) 
Substitute (J23J),(|2U) and (EH in (DUD putting m = 0. 

/So r n-1 ^ 

2/o(^) = c ^m^u + i U 2n+A (26) 

n=0 Ui=0 J 

Observe the terms inside curly brackets of sequence c n which include one term of A n 's in 
( 112"]) : c n with odd index (c\, C3, 05,- • ■ ). 



(a) As /3 =0, then B 1= in (1T2]) . 



c i — A) c o 

c 3 = Ao-S 2 c 
c 5 = A B 2 ,4C 

(27) 

C7 = A B2,4fiC 
C9 = ^0-^2,4,6,80) 
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As i=l in (1221 . 



B Wl+2 = where ft = 0, 1, 2, • ■ • (28) 
Substitute (|27jl in ( ITU]) putting m = 1 by using ([28]) . 



r«-i 

!/?(*) = c A £ J] B 2h+2 \ * 2 "'+ 1+A (29) 



n=0 Ui=0 



In fl29|) y°(2;) is sub-power series, having sequences c n including one term of A n 's in ([7]) as 
Po=0, for the polynomial case in which makes B n term terminated. 

(b) As /? =1, then B 3 =0 in (fT2"l). 



Ci — ^4o c o 

C3 = {A 5 2 -l + A 2 -lfi 1 }co 

c 5 = {A 5 2i4 • 1 + A 2 5 4 fii}c 
c 7 = {A 5 2i4 ,6 • 1 + A 2 B 4i6 B 1 }c Q 
c 9 = {A B 2Afii8 • 1 + A 2 B 4fij8 B 1 }c 



(30) 



The first term in curly brackets of sequence c n in (130]) is same as (1271) . Then, its solution is 
equal to (129]) . Substitute (128]) into the second term in curly brackets of sequence c n in (130|) . 
Its power series expansion including the first and second terms in curly brackets of sequence 
c n in (13"0"|) . analytic at x — 0, is 



ft r«-i ^ ft r«-i 

yiCn) = co \A J2 \U B 2n+2 + A 2 B 1 £ \ ]J B 2ll+2 \ \ x 2n+l+x (31) 

n=0 Ui=0 J n=l Ui=l 



In (13T|) is sub-power series, having sequences c n including one term of A n 's in ((TJ, as 

/?o=l for the polynomial case in which makes B n term terminated. 
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(c) As (3 =2, then B 5 =0 in (JTJ). 



ci = A c 

c 3 = {A B 2 ■ 1 + A 2 ■ 1Bi}cq 

c 5 = {A B 2A ■ 1 + A 2 B 4 B 1 + A 4 ■ 15 lj3 }c 

C 7 = {A)5 2 ,4,6 • 1 + ^B^B, + A 4 B 6 B^}c (32) 
C9 = {A)^2,4,6,8 • 1 + A 2 B Afi:8 B 1 + A 4 B 6t8 B 1 ^}c 
en = {^-0-82,4,6,8 ■ 1 + ^-2-84,6,8-81 + A 4 I? 6i8) io-Bi,3}co 

The first and second term in curly brackets of sequence c n in (1321) is same as (I30I) . Then its 
power series expansion is same as fl3Tj) . Substitute (128]) into the third term in curly brackets 
of sequence c n in (1321) . Its power series expansion including the first, second and third terms 
in curly brackets of sequence c n in (I3"2"j) . analytic at x — 0, is 



2n+l+A 



{ft f n-1 ^ ft r n-1 ^ ft ( n-1 ~\ 

A °J2\ l[B2 il+ A+A 2 B 1 J2\ llB^+A+AiB^J^i 11^ 
n=0 I ii=0 J n=l I u = l J n=2 I i 1= 2 ) ) 

(33) 

In (133]) (x) is sub-power series, having sequence c n including one term of A n 's in ([7]) as 
Po=2, for the polynomial case in which makes B n term terminated. By repeating this process 
for all fto — 3, 4, 5, • • • , I obtain every y{(x) terms where j > 3. According to f[2"Uj) . f[3"Tj) . f l3"3"j) 
and every t/j(x) where j > 3, the general expression of yi(x) for all (3q is 



ft r »o-i ft r*2-i 

yi (x) = co j2\ n e n [ i> ^ 2j2+i+a (34) 

io=0 L ii=0 «2=io L«3=*0 



Observe the terms of sequence c n which include two terms of A n 's in ([15]) : c n with even 
index (c 2 , c 4 , c 6 ,- ■ ■ ). 
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(a) As (3 =0, then B 1= in 05). 



C2 = A ,ic 

c 4 = A {A! • 1 • B 3 ■ 1 + A 3 B 2 ■ 1 ■ l}c 

c 6 = A {A! • 1 • B 3 , 5 • 1 + A 3 B 2 B 5 ■ 1 + A 5 B 2A ■ 1 • l}c 

c 8 = A {A 1 ■ 1 ■ 53,5,7 ■ 1 + A 3 B 2 B 5J ■ 1 + A 5 B 2A B 7 ■ 1 + A 7 B 2A& ■ 1 ■ l}c 

Cio = ■ 1 " £3,5,7,9 " 1 + ^3-82-85,7,9 • 1 + ^5-82,4-87,9 ■ 1 + A 7 B 2 ^ 6 B 9 ■ 1 + y4 9 5 2 ,4,6,8 ■ 1 " l}co 

(35) 

(i) As /3i=0, then B 2 =0 in (J35]). 

C2 = v4 ,lC 

c 4 = v4 Ai • 1 • B 3 - lc 
c 6 = A A 1 ■ 1 • £3.5 • lc 

(36) 

C 8 = AqAx ■ 1 • £3,5,7 • lc 
Cio = AqAi ■ 1 • -B 3 ,5,7,9 • lCo 

As i=2 in (1221), 

5 2/32+ 3 = where & = 0,1,2,- •• (37) 
Substitute ( 13"B"1) in ( 1TU|) putting m = 2 by using (IBTj) . 

& r n-1 ^ 

y°/{x) = c A A 1 II ^u+3 * 2n+2+A (38) 

n=0 Ui=0 J 

In (138]) y^' ^) 1S sub-power series, having sequences c n including two term of A^s in (JTj) as 
/3o=0 and /?i=0, for the polynomial case in which makes B n term terminated. 
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(ii) As J3x=l, then 54=0 in flM). 

C 2 = A 0A c 

c 4 = AoiA, • 1 • # 3 • 1 + A 3 # 2 • 1 • l}c 
c 6 = A {A! • 1 • £3,5 ■ 1 + A 3 B 2 B 5 ■ l}c 

(39) 

c 8 = AoiA, ■ 1 ■ £3,5,7 ■ 1 + A 3 B 2 B 5J ■ l}c 
cio = A){^i ■ 1 ■ #3,5,7,9 • 1 + A 3 B 2 B 5Jt9 ■ l}c 

The first term in curly brackets of sequence c n in ( l39l) is same as ( 1361) . Then its power 
series expansion is equal to (1381) . Substitute ( 1371) into the second term in curly brackets of 
sequence c n in ( 139]) . Its power series expansion including the first and second terms in curly 
brackets of sequence c n , analytic at x = 0, is 

{fa f n-l ^ fc fn-1 ^ 

a ^ e n b ^ + ^ ^ n b ^ x2n+2+x ( 4 °) 
n=0 Ui=0 J n=l Ui=l J J 

In (HOI 2/2' C^) i s sub-power series, having sequences c n including two term of A n 's in ((7j) as 
/3o=0 and /?i=l, for the polynomial case in which makes B n term terminated. 

(iii) As /3i=2, then # 6 =0 in 

C2 = A),1 C 

C 4 = A){^1 ■ 1 ■ S 3 ■ 1 + A 3 # 2 ■ 1 • l}c 

C 6 = A {A 1 • 1 • S 3 ,5 • 1 + A 3 B 2 B 5 • 1 + ^5^2,4 • 1 • l}c 

(41) 

c 8 = A {A 1 ■ 1 ■ #3,5,7 • 1 + A 3 B 2 B 5J ■ 1 + A 5 B 2A B 7 ■ l}c 

Cio = A {A 1 ■ 1 • #3,5,7,9 • 1 + A 3 # 2 #5,7,9 " 1 + A 5 #2,4#7,9 " l}c 

The first and second term in curly brackets of sequence c n in (I4T1) is same as (|39l) . Then its 
power series expansion is same as ( 1401) . Substitute ( 137|) into the third term in curly brackets 
of sequence c n in ( 14"TT) . Its power series expansion including the first, second and third terms 
in curly brackets of sequence c n in ( |4T|) . analytic at x = 0, is 

r #2 r n-l ] & ( n-l ^ P a ( n-l ^ ^ 

2/2 o ' 2 (x)=c Ao^ 1 ^i n^ 1+ 3 +^2^^ n ^+3 +^2,4 n^ 1+ 3 U 2n+2+A 

I n=0 I ii=0 J n=l I ii=l J n=2 I ii=2 J J 

(42) 
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In (T4"2j) y 2 ,2 (x) is sub-power series, having sequences c n including two term of A n 's in (j7j) as 
/3o=0 and /3i=2, for the polynomial case in which makes B n term terminated. By repeating 
this process for all (3\ = 3, 4, 5, • • • , we obtain every y 2 '\x) terms where j > 3. According to 
fl38|) . fj40|) . f )42|) and every y 2 ,:i (x) where j > 3, the general expression of y 2 (x) for the case 
of (5q = replacing the index n by i 9 is 

01 ( io-l 02 ( %2 i ^ ^ 

y 2 °(x) = c A A 2i0+1 B 2ll+2 £ II ^ +2+A (43) 

io=0 L ii=0 i2=io Ki3=io ) ) 



In fj4"3"j) ( x ) * s sub-power series, having sequences c n including two term of A n 's in ((7]) as 
/?o=0, for the polynomial case in which makes B n term terminated. 

(b) As A,=l, then S 3 =0 in (JTHD - 



C 2 = AqiCq 



c 4 = j A 
c 6 = {A) 
c 8 = {A 
+A 2 
cio = <^ A 



A x ■ 1 • 5 3 • 1 + A 3j B 2 -1-1 



+ A 2 



A 3 ■ 1 ■ 1 ■ Si 



Co 



A x • 1 • S3, • 1 + A 3 B 2 B 5 ■ 1 + A b B 2A -1-1 



+ A 2 



As-l-BsBi+AsB^l-B^co 



Ax • 1 • £3,5.7 • 1 + A 3 B 2 B, 7 ■ 1 + A 5 B 2 aB 7 ■ 1 + A 7 B 2 , 6 -1-1 



A 3 • 1 • B, 7 B 1 + AsB±B 7 B x + A 7 B, 6 • 1 ■ Bx 



Cq 



Ai • 1 • -63,5,7,9 " 1 + A 3 S 2 55,7,9 • 1 + A 5 5 2 ,4i?7,9 ■ 1 + A 7 B 2A fiB 9 ■ 1 + A 9 i? 2 ,4, 6 , 8 ■ 1 ■ 1 

+A 2 A 3 ■ 1 ■ #5,7,9!?! + A 5 B i B 7t9 B 1 + A 7 B ifi B 9 Bi + A 9 S 4 , 6 , 8 ■ 1 ■ i?x |c 



(44) 

The first square brackets including Ao inside curly brackets in sequence c n in (144"1) is same as 
( |35l) . Then its power series expansion is same as ( 1431) . Observe the second square brackets 
including A 2 inside curly brackets in sequence c n in (j44 



(i) As then £? 4 =0 in the second square brackets including A 2 inside curly brackets 
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m sequence c„ m ^ , 

c 4 = A 2 B 1 {a 3 • 1 • l}c 

c 6 = A 2 B^A 3 ■ 1 ■ B 5 y o 

c 8 = A 2 S 1 {a 3 -1- j B 5i7 }c (45) 

Cio = A 2 B 1 ^A 3 ■ 1 • #5,7,9 j C 

Its power series expansion of flU]) by using (l"T"*j) . analytic at x = 0, is 

& r n-1 ^ 

^(x) = c A 2 B 1 A 3 II ^ 1+ 3 x M+A (46) 

n=l Ui=l J 

In ( 146]) T/g' 1 ^) is sub-power series, for the second square brackets inside curly brackets in 
sequence c n including two term of A n 's in flUJ) as A)=l and /?i=l, for the polynomial case 
in which makes B n term terminated. 

(ii) As /3i=2, then B e =0 in second square brackets inside curly brackets in sequence c n 
including A 2 in (jUJ. 

C4 = v4 2j B 1 |a 3 -1-i}co 

c 6 = AsSxjAg • 1 • B 5 + A 5 B 4 ■ l}c 

c 8 = A 2 B 1 {A 3 ■ 1 ■ #5,7 + A 5 B 4 B 7 y o (47) 



cio = A 2 Bx{a 3 ■ 1 • 55,7,9 + A 5 B 4 B 7 ^c 



The first term in curly brackets of sequence c n in (147"]) is same as fj45l) . Then its power series 
expansion is same as (14*6]) . Substitute (I37"j) into the second term in curly brackets of sequence 
c n in ([47]) . Its power series expansion including the first and second terms in curly brackets 
of sequence c n in (147j) . analytic at x = 0, is 



/?2 f n.— 1 ^ fii ( n—1 



y l 2 \x) = c Q A 2 B x \A 3 J2 \ II 5 ^+ 3 f + ^ E \ II 

^ n =l Ui=l J ri=2 ln=2 

,1,2/ 



-3 > i> X 2 " +2+A 



In (145]) y 2 ' (x) is sub-power series, for the second square brackets inside curly brackets in 
sequence c n including two term of A n 's in (1441) as /?o=l and /3i=2, for the polynomial case 
in which makes B n term terminated. 
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By using similar process as I did before, the solution for /?i=3 with Bg=0 for the second 
square brackets inside curly brackets in sequence c n including A 2 in (TJ4]) is 



02 (n 1 ^ 02 fn-l ~\ 02 (n 1 

n=l l.ii = l ) n=2 l_ii=2 ) ra=3 l_ii=3 



+3 



(49) 

By repeating this process for all (5\ = 4, 5, 6, • ■ ■ , we obtain every y\ 3 (x) terms where j > 4 
for the second square brackets inside curly brackets in sequence c n including two term of 
A n 's in (T441) . According to (T43]), fj46l) . (JHJ), fT49|) and every y 2 J (x) where j > 4, the general 
expression of y\ (x) for all /3 = 1 replacing the index n by i is 



0i ( *o-l /8a r«2-l 

j/£(x) = c <{ A S A 2i +i 11 B2i ^+ 2 \ n 

io=0 L ii=0 i2=«o L«3=io 

01 ( io — 1 fti f 12 — 1 



+^1^ ]T J A 2i0+1 H B 2il+2 J2\Ub 



2i 3 +3 



io=l 



22=«0 V«3=*0 



(50) 



>x 



2%2+2+X 



(c) As /3 =2, then ,85=0 in (fl5]) . 



C2 — A),i c o 



C 4 = |A • 1 
C 6 = {A • 1 

+AiB h3 



A 1 -l-B 3 + A 3 B 2 ■ 1 



A 2 B X 



A, ■ 1 ■ 1 



Ax • 1 • S 3i5 + A 3j B 2j B 5 + A 5j B 2 , 4 • 1 



+ A 2 B X 



A 3 ■ 1 • B 5 + A b B A ■ 1 



As ■ 1 • 1 



Co 



r-.s = |.l n • 1 

cio = |A • 1 

+a 4j b 1i3 



Al • 1 • -83,5,7 + A 3 .B 2 £?5 )7 • 1 + A 5 _B 214 i? 7 + A 7 -B 2j 4 )6 • 1 



A 3 • 1 ■ 5 5j7 + A 5 BiB 7 + A 7 B ifi ■ 1 



+ ; 



A 5 " 1 ■ E 7 + A 7 5 6 ■ 1 



c 



Ai • 1 • -83,5,7,9 + A 3 B 2 B 5 j$ + A 5 B 2ti B 7j9 + A 7 B 2ti fiB 9 + A 9 i? 2i 4 )6i 8 • 1 
A3 • 1 ■ -85,7,9 + A5_B45 7i 9 + A 7 £? 4)6 I?9 + A 9 -B 4i6i8 ■ 1 



A5 • 1 • B 7 9 + A 7 B 6 B g + A 9 S 6 g • 1 



Co 



(51) 

By repeating similar process from the above, the general expression of y 2 (x) for all (3q = 2 
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in (EH) is 



{ft ( io-l ft 

a o e \ A2i °+ x n b ^+ 2 e \ n b ^ +3 
i =0 L ii=0 «2=«0 l«3=*0 

ft f i -l ft fi 2 -l 

+A 2 B X £ \ A2i 0+ i n S ^ 1+ 2 £ II S ^3 + 3 }> }> (52) 

i =l L U=l i2=«0 L«3=*0 

ft r i -i ft r i 2 -i 

+A 4 Si i3 £ S ^2i +l Y\_ B 2h+2 £ S H ^3+3 
«0=2 I il=2 «2=«0 V*3=»0 



>x 2j 2 +2+A 



Again by repeating this process for all (3q = 3,4,5, ••■ , I obtain every y J 2 {x) terms where 
j > 3. Then I have general expression y2{x) for all (3q of two term of A n 's according fT43|) . 
( 150|) . (|52|) and 2/2 ( x ) terms where j > 3. 

A> r io-i ft r »2-i ft r? 4 1 ^ ^ ^ 

^(x) = c £ ^ ^ n e ^ 2+ i n ^3+2 £ ^ n ^ 5+3 * 2i4+2+A 

io=0 I. ii=0 «2=io I «3=*0 14=12 1*5=22 ) ) ) 

(53) 

By using similar process for the previous cases of zero, one and two term of A ra 's, the function 
ys(x) for the case of three term of A^s is 

i — 1 ft f i 2 -l 

y 3 (z) = c > J <( A 2io Yl B 2ll +i £ \ A 2i 2 +i II B ^+ 2 ( 54 ) 

u=0 «2=«o I h=io 

t4-l ft fi 6 -l 1111 

> n s - 5+3 £ n ^ 7+4 k 6+3+A 

«5=«2 «6=«4 L*7=*4 ) ) ) ) 

By repeating this process for all higher terms of A^s, I obtain every y m {x) terms where 
m > 3. Substitute ( 126|) . ( 13~4"|) . (|53|) . (154"]) and including all y m (x) terms where m > 3 into (|9]). 
The general expression of y(x) for the polynomial case which makes B n term terminated in 
the three term recurrence relation is 

{ft /io-l \ ft ( io-l ft /i 2 -l \ 

e n ^° +a + e 4* n £ n ^ 3+2 ^ 2+i+a 

io=0 \ii=0 / io=0 L ii=0 «2=«o \i3=*0 / J 

00 ( ft r io-l iV-1 / ft i 2fe -l 

+ e s e 1 n s 2n+i n ( e ^w+fc 11 ^2i 2fc+1 +(fc+i) 

V=2 I i =0 L ii=0 k=l \«2fe=«2(fe-i) «2fc+l=»2(fc-i) 

x E ( II ^ +1+ (^d) ^^ +iV+A ] (55) 



ft 1 




E 


^2i 


i =0 I 




ft 




£< 

14=12 


|^2i 4 



«2JV=«2(AT-1) \ *2iV+l=«2(jV-l) 
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For infinite series, replacing j3o,(3i,/3k and (3^ by oo in (1551) 

{oo /*0— 1 \ oo (■ i — 1 oo /i2— 1 \ 

^ n b» 1+1 x 2j ° +A + £ | a 210 n 5 2il+1 ^ n s 2l3+2 j> x 
i =0 Vi=Q / io=0 I ii=0 «2=io V3=*0 / J 

OO f OO f IQ — 1 TV— 1 / OO »2fc — 1 

JV=2 L i =0 I h=0 k=l \»2fc=*2(fc-i) *2fc+i=«2(fc-i) 



2i 2 +l+A 



OO / «2JV — 1 

X 



e n ^ (56) 

*2JV=i2(iV-l) V *2JV+1=«2(JV-1) / ) ) ) 

is exactly equivalent to ([2"T]) . (151)]) is the another general expression of for the 
infinite series. 

V. RECURRENCE RELATION AND FIBONACCI NUMBERS OF 
HIGHER ORDER 

A. Three term recurrence relation and Fibonacci sequence 

The Fibonacci sequence is: 

1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, • • • (57) 
The numbers in the sequence follows the recursive relation 

c„+i = c n + c n _i : n > 1 (58) 

with seed values 

c = 1 ci = 1 (59) 
The power series of the generating function of the Fibonacci sequence is 



x 

c n x - — 

n=0 

As we know, the three term recurrence relation is 



X>x"= X - (60) 
^— ' 1 — a; — ar 



c n+ i = A n c„ + B n c n _i ; n > 1 (61) 

with seed values 

ci = A c (62) 
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© is the expansion of (I6TI) . If coefficients A n = B n = 1 in (16T|) and Co = ci = 1 in (|62l) . 
then c n follows Fibonacci sequence. Lucas series is another sequence generated by the three 
term recurrence relation with constant coefficients A n = B n = 1 in ( 16 ip and Cq = 2,c\ = 1 
in (16"2"|) . You can think of these two sequences having constant coefficients as the most basic 
three term recurrence relation. In contrast, Heun and Mathieu equations coefficients A n and 
B n are defined to be non-constant (second order polynomial in denominator and numerator: 



see (3), (4a)-(4c) in Ref. [13| and (4), (5a)-(5c) in Ref. |15|) The generating function of the 
Fibonacci sequence corresponds to infinite series of three term recurrence relation. If all 
coefficients A n = B n = 1 in f[2TT) and (151)1) . then it is equivalent to ( 16TJ]) . 

B. Two term recurrence relation 

There is an algebraic number sequence in which is 

1,1,1,1,1,1,- •• (63) 
I call (163]) the identity sequence, and it's recurrence relation is 

c n +i = c n : n > (64) 

with seed values 

c = 1 (65) 
The power series of the generating function of the identity sequence is 

oo 1 

$> n x" = - (66) 

If B n =0 in flBTj) . then three term recurrence relation becomes a two term recurrence relation. 

c, t+1 = A n c n ; n > (67) 

Some examples are the Legendre function, Kummer function, hypergeometric function, 
Bessel function, etc. The number of each sequence c n in ( 167|) is 

1,1,1,1,1,1,- •• (68) 

( 168]) is equivalent to ( 163]) . As I put A n =l in (l6"Tj) . it becomes (1641) . You can think of this 
sequence having constant coefficients as the most basic two term recurrence relation. The 
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power series expansion of ( |67|) for the infinite series is 



oo oo /n—1 \ 

y(x) = Cnx n+X = co II ^ ^ ( 69 ) 

n=0 n=0 \i=0 / 

And polynomial case of f l6T|) is 

y(s) = ^ c n a;" +A = Co £ J] A, x" +A (70) 

n=0 n=0 \i=0 / 

The generating function of the identity sequence corresponds to infinite series of two term 
recurrence relation. If all coefficients A n = 1 in ( 1691) . then it definitely equivalent to fl6"6"|) . 

C. Four term recurrence relation and Tribonacci sequence 

Now, let's think about four term recurrence relation in ordinary differential equation. 
The four term recurrence formula is 

c n+l = A n C n + B n C n _i 

+ C n c n _ 2 ; n > 2 (71) 

with seed values 

ci = A c c 2 = (A A 1 + Si) c (72) 
And the number of each of sequence c n in (I7TT) is the following way: 

1,1,2,4,7,13,24,44,- •• (73) 

(173p is Tribonacci number, and it's recurrence relation is 

c„+i = c n + c„_i + c„_ 2 ; n > 2 (74) 

with seed values 

c = d = 1 c 2 = 1 (75) 

If A n = B n = C n = 1 in (j71j) . it's exactly equivalent to Tribonacci recurrence relation. Four 
term recurrence relation (non-constant coefficients A n , B n and C n ) is the more general form 
than Tribonacci recurrence relation (constant coefficients A n , B n and C n ). 
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D. Five term recurrence relation and Tetranacci sequence 



And five term recurrence relation in ordinary differential equation is 

Cn+1 An Cn ~t~ Cn—1 + CnC n _ 2 + -D n c n _ 3 ; n > 3 (76) 

with seed values 

Cl = A co, c 2 = (A A! + B x ) c , c 3 = (AoAiAs + A 5 2 + A 2 5i + C 2 ) c (77) 
And the number of each of sequence c n in (1761) is the following way: 

1,1,2, 4, 8, 15, 29, 56, 108, 208, ■ • • (78) 
( ITS"]) is Tetranacci number, and it's recurrence relation is 

c n +i — c„ + c n _i + c„_2 + c n _3 ; n > 3 (79) 

with seed values 

c = ci = 1 c 2 = 1 c 3 = 2 (80) 

If A„ = £> n = C„ = D n = 1 in (176|) . it becomes exactly equivalent to Tetranacci recurrence 
relation. From the above multi-term recurrence relation is more general form than n-nacci 
recurrence relation. 

VI. CONCLUSION 

In this paper I show how to generalize three-term recurrence relation for polynomials 
and infinite series analytically. In the next papers I will work out the analytic solution 
for the three term recurrence relation such as Heun equation and its confluent form, Lame 
equation, Mathieu equation. I will derive the power series expansion, the integral form and 
the generating function of Heun, Mathieu, etc functions, (see section VII for more details) 

VII. SERIES "SPECIAL FUNCTIONS AND THREE-TERM 
RECURRENCE FORMULA (3TRF)" 

This paper is 2nd out of 10. 

1. "Approximative solution of the spin free Hamiltonian involving only scalar potential 
for the q — q system"— - In order to solve the spin-free Hamiltonian with light quark masses 
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we are led to develop a totally new kind of special function theory in mathematics that 
generalize all existing theories of confluent hypergeometric types. We call it the Grand 
Confluent Hypergeometric Function. Our new solution produces previously unknown extra 
hidden quantum numbers relevant for description of supersymmetry and for generating new 
mass formulas. 

2. "Generalization of the three-term recurrence formula and its applications" 12 - Gener- 
alize three term recurrence relation in linear differential equation. Obtain the exact solution 
of the three term recurrence relation for polynomials and infinite series. 

3. "The analytic solution for the power series expansion of Heun function" 13 - Apply 
three term recurrence formula to the power series expansion in closed forms of Heun function 
(infinite series and polynomials) including all higher terms of A n s. 

4. "Asymptotic behavior of Heun function and its integral formalism" ,— - Apply three 
term recurrence formula, derive the integral formalism, and analyze the asymptotic behavior 
of Heun function (including all higher terms of A n s). 

5. "The power series expansion of Mathieu function and its integral formalism" , 15 - Apply 
three term recurrence formula, analyze the power series expansion of Mathieu function and 
its integral forms. 

6. "Lame equation in the algebraic form"-i£ - Applying three term recurrence formula, 
analyze the power series expansion of Lame function in the algebraic form and its integral 
forms. 

7. "Power series and integral forms of Lame equation in the Weierstrass's form and its 
asymptotic behaviors"— - Applying three term recurrence formula, derive the power series 
expansion of Lame function in the Weierstrass's form and its integral forms. 

8. "The generating functions of Lame equation in the Weierstrass's form"— - Derive the 
generating functions of Lame function in the Weierstrass's form (including all higher terms 
of A n 's). Apply integral forms of Lame functions in the Weierstrass's form. 

9. "Analytic solution for grand confluent hypergeometric function"— - Apply three term 
recurrence formula, and formulate the exact analytic solution of grand confluent hyperge- 
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ometric function (including all higher terms of A n 's). Replacing fi and ecu by 1 and — q, 
transforms the grand confluent hypergeometric function into Biconfluent Heun function. 

10. "The integral formalism and the generating function of grand confluent hypergeomet- 
ric function"— - Apply three term recurrence formula, and construct an integral formalism 
and a generating function of grand confluent hypergeometric function (including all higher 
terms of A n 's). 
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